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Abstract. This article is dedicated to the memory of Vadim Kuznetsov, and begins with 
some of the author's recollections of him. Thereafter, a brief review of Somos sequences is 
provided, with particular focus being made on the integrable structure of Somos-4 recur- 
rences, and on the Laurent property. Subsequently a family of fourth-order recurrences that 
share the Laurent property are considered, which are equivalent to Poisson maps in four 
dimensions. Two of these maps turn out to be superintegrable, and their iteration furnishes 
infinitely many solutions of some associated quartic Diophantine equations. 
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1 Introduction 

It is with considerable sadness that I begin to write this piece in memory of Vadim Kuznetsov, 
whose death came as a great shock to me. However, I do not wish to remain in melancholic 
mode, but rather I would like to recall some of my fondest and happiest memories of him. 

While I was a PhD student in Edinburgh, I used to travel to Leeds every so often to attend 
the LMS workshops on integrable systems, and I'm sure I must have first met Vadim at one 
of these meetings. To begin with I remember his charming smile, as well as his relaxed way of 
asking penetrating mathematical questions. I also recall the great enthusiasm and energy with 
which he would give a seminar, and his clarity of presentation. 

Shortly after graduating from Edinburgh, in September 1997 I went to Rome to take up my 
first postdoctoral position, working with Orlando Ragnisco in the Physics Department of Roma 
Tre. It was during this period that I had the privilege of getting to know Vadim a lot better. 
A few months after my arrival, he came to Rome to visit Orlando for a month, and the three of 
us ended up working together on a project that was suggested by Vadim, concerning Backlund 
transformations (BTs) for finite-dimensional integrable Hamiltonian systems. This turned out 
to be very fruitful, resulting in three joint publications [231 EH E^- 

Vadim's presence in Rome was immensely stimulating for me, because he succeeded in posing 
just the right question, at the precise moment when I had the necessary tools available to 
answer it. The specific problem that he first presented to me and Orlando was the construction 
of BTs for certain integrable classical mechanical systems corresponding to reduced Gaudin 
magnets. A particular concrete example of such a system was the case (ii) Henon-Heiles system, 
an integrable system with two degrees of freedom. As it happened, in my PhD thesis I had 
already constructed an analogous BT for the non-autonomous case of this system, as well as 
deriving the explicit formula for the generating function of the canonical (contact) transformation 
in that case [22]. During my viva voce examination a few months earlier, Allan Fordy had 
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actually asked me whether the same sort of derivation could also be applied to the autonomous 
case, to produce a Poisson correspondence in the spirit of [TB], and I could see no obvious 
obstruction. Thus it was that, when Vadim arrived in Rome, his vivid explanation of BTs, as 
well as his insistence that we should start constructing new ones, was all that I needed to work 
out the BT for the Henon-Heiles system |23] , and this soon revealed a similar algebraic structure 
underlying many other examples |24j . 

After he left Rome, I saw Vadim again in June 1998 at the conference Integrable Systems: 
Solutions and Transformations in Guardamar, Spain, where he came with his wife, Olga, and 
his son, Simon. We sat down together in the sunshine and completed some of the work on 
the second paper [23] while we were there. Subsequently, I saw Vadim sporadically at various 
meetings in Leeds and elsewhere, and we always found the time for a friendly chat about our 
lives and work. I particularly remember a very brief and enjoyable (but fiery) dispute that 
we had in Cambridge in 2001, while sitting together during an interlude between lectures in 
the Newton Institute. It boiled down to a minor difference in our points of view, which we 
respectively argued for without compromise, so that (having each seen the other's perspective) 
there was no love lost between us. 

The rest of this article is concerned with a family of discrete dynamical systems (Poisson 
maps) in four dimensions, the first few of which are integrable, while the others are not. Before 
going into details, I should like to explain why I have chosen this topic. The work I did in my PhD 
was primarily concerned with integrable systems in the continuous setting (ordinary and partial 
differential equations), and it was not until Vadim's visit to Rome that I began to get actively 
interested in discrete systems. Ever since then, I have found the subject of discrete dynamics 
increasingly fascinating, and I shall always have Vadim to thank for inspiring me to look in this 
direction. Another interesting and unexpected property of the Poisson maps considered below is 
that their iterates are Laurent polynomials in the initial data; this is an instance of the Laurent 
phenomenon [T3]. Vadim was an expert on special functions, and orthogonal polynomials in 
particular (for one of his many contributions in this area, see [35], for instance). However, most 
of the sequences of (Laurent) polynomials treated below satisfy nonlinear equations instead of 
linear ones. 

The theory of discrete integrable maps has seen a great deal of activity in the past twenty 
years. The situation was much clarified by Veselov [571 EH [59] who introduced integrable 
Lagrange correspondences - a natural discrete-time analogue of Liouville integrable continu- 
ous flows - which induce (generically multi-valued) shifts on the associated Liouville tori (see 
also [3]). Given a continuous integrable system, it is natural to seek a discretization of it that 
retains both the integrability and as many other properties as possible (e.g. Poisson structure, 
Lax pair, etc.). However, in general such a time-discretization will be implicit, and it will not 
preserve the same integrals as the original continuous system (see [52] for the state of the art 
in integrable discretizations). Building on results obtained for the Toda lattice by Pasquier 
and Gaudin [38], Kuznetsov and Sklyanin identified a special class of time-discretizations for 
integrable Hamiltonian systems which they referred to as BTs [32], by analogy with Backlund 
transformations for evolutionary PDEs. 

In the setting of finite-dimensional systems with a Lax pair, BTs were identified as explicit 
Poisson maps which preserve the same set of integrals as the continuous flow that they discretize, 
and depend on a Backlund parameter A which satisfies a certain 'spectrality' property with 
respect to a conjugate variable (where (A, /j.) are the coordinates of a point on the spectral 
curve associated with the Lax pair). The viewpoint that I emphasized in [23^ [24} [26] was 
that the systems being considered were reduced/stationary flows of the KdV hierarchy, whose 
BTs could be obtained by reduction from the Darboux-Backlund transformation for KdV (this 
is in the same vein as the dressing chain [60] - see also |63j). while the BTs in [2^ were 
derived more directly. In Vadim's work with Pol Vanhaecke [34j . all of the previously known 
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examples were unified via an algebro-geometric approach, which explained the deeper meaning 
of BTs as discrete shifts on the (generalized) Jacobian of the associated spectral curve, thus 
identifying them as the discrete-time counterparts of algebraically completely integrable systems, 
as described in [56], for instance. While there has been subsequent work by Vadim and others 
on BTs in classical mechanics [H \T3\ [39] , a lot of the original motivation for studying them 
came from quantum integrable systems (Baxter's Q-operator). This idea has proved extremely 
effective (see e.g. [MIES]), and will no doubt continue to bear fruit for a long time to come. 

The last time I saw Vadim was in Leeds in April 2005, when he invited me to give one 
in the series of Quantum Computational seminars that he organized ther^. At that time 
I spoke about Somos sequences, which are reviewed in the next section. Li the evening after 
the seminar I went out for a very enjoyable dinner with Vadim and Olga, together with Oleg 
Chalykh and Sara Lombardo. I made an appointment to see Vadim in his office early the next 
morning, so that before my return home we had a good discussion about his recent work on the 
integrable dynamics of spin chains that arise in models of Fermi-Bose condensates [65] and BCS 
superconductors [66], and he described an unsolved problem concerning special solutions. This 
is how I remember him now: full of energy and always seeking to answer new questions. 



2 Somos sequences and the Laurent property 

The properties of integer sequences generated by linear recurrences have been the subject of 
a great deal of study in number theory, and nowadays they find applications in computer science 
and cryptography [12]. However, the theory of nonlinear recurrence sequences is still in its 
infancy. Clearly, a fcth-order nonlinear recurrence relation of the form 

^n+k ~ FiXfi^ Xn-\-1i ■ ■ ■ l) (2-1) 

is just a particular sort of discrete dynamical system, so such recurrences can be considered 
as generating a special type of nonlinear dynamics. If we want (j2.1|) to generate sequences of 
integers, then choosing F to be a polynomial with integer coefficients will certainly do the trick, 
but in general the corresponding map in R'"' (or C*^) will not have a unique inverse. Moreover, 
in that case such sequences generically exhibit double exponential growth i.e. loglx^l grows 
exponentially with n. A simple example in this class is the quadratic map defined by the 
recurrence 

Xn+i = xl + c (2.2) 

with a parameter c, which is a prototypical model of chaos. However, note that the special 
cases c = 0,-2 are exactly solvable [9], and in these cases one can also argue that ()2.2p is 
integrable in the sense of admitting a commuting map (see [57] and references). The theory of 
linear recurrence sequences relies heavily on the fact that they are explicitly solvable. Thus it 
is natural to look for nonlinear recurrences that share this property, or that are integrable in 
a broader sense. 

In the case that the map corresponding to ()2.ip is invertible, one can also allow F to be 
a rational function, thereby considering birational maps, but then it is no longer clear that 
integer sequences should result. However, it turns out that among those rational recurrences of 
the particular form 

Xn+k — f{Xn+l-, ■ • • ■> Xn+k—l)i (^-3) 
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there is a very large class of recurrences that generate integer sequences from suitable initial 
data. One of the first known examples of this type is the Somos-4 recurrence 

Xn+A Xn = a Xn+Z Xn+1 + P {Xn+2f , (2.4) 

which was found by Michael Somos when he was investigating the combinatorics of elliptic theta 
functions. Somos observed numerically that by taking the coefficients a = (3 = 1 and initial 
data xo = xi = X2 = X3 = 1, the fourth-order recurrence ()2.4p yields a sequence of integers (50] . 
that is 

1, 1, 1, 1, 2, 3, 7, 23, 59, 314, 1529, 8209, 83313, .... (2.5) 
Similarly he noticed that for the Somos-/c recurrences 

Xn+k Xn — ^ '] CXj Xn+k—j Xn+j (2 '6) 

i=i 

with all coefficients aj = 1, if all k initial values are 1 then an integer sequence results for 
k = 5,6, 7, but denominators appear for A: = 8. 

Various direct proofs that the terms of the sequence (|2.5p are all integers were found at 
the beginning of the 1990s, when various other examples were found [T71 STj, but a deeper 
understanding came from the realization that the recurrence (j2.4p has the Laurent property: its 
iterates are all Laurent polynomials in the initial data (and in a, (3) with integer coefficients. To 
be more precise, the iterates of (j2.4p satisfy Xn G l^lx^^ ,x^^ ,X2^ ,x^^ ,a, (3] for all n, from which 
the integrality of the particular sequence (j2.5p follows immediately. A little earlier, when Mills, 
Robbins and Rumsey made their study of the Dodgson condensation method for computing 
determinants |36j (which produced the famous alternating sign matrix conjecture [2]), they 
considered the recurrence 

De^rn,n+lDe^rn.^n-~l = « De+i^rn^nDi-i^rn,n + P De^rn+l,nDe^rn-l,n, (2.7) 

for a = 1 and observed that it produced Laurent polynomials in the initial data. The equa- 
tion (j2.7p thus became known within the algebraic combinatorics community, where it is referred 
to as the octahedron recurrence |43j . while in the theory of integrable systems it is known 
as a particular form of the discrete Hirota equation [68] (the bilinear equation for the tau- 
function of discrete KP). The Somos-4 recurrence (12. 4p is an ordinary difference reduction of 
the partial difference equation (j2.7p : it has been noted by Propp that if Xn satisfies (j2.4p then 
De^m,n = X2n+m Satisfies the discrete Hirota equation (see also [51] for another reduction). 

Many more examples of this Laurent property have begun to emerge quite recently as an 
offshoot of the theory of cluster algebras due to Fomin and Zelevinsky (see [15] and references) . 
The exchange relations in a cluster algebra of rank k are typified by a recurrence of the form 

Xn+kXn = Cl Mi(x„+i, . . . ,Xn+k-l) + C2 JVhiXn+l, ■ ■ ■ , Xn+k~l) (2.8) 

for suitable monomials Mj and coefficients Cj, which is a special case of (j2.3p . In [14j . the general 
machinery of cluster algebras was shown to be very effective in proving the Laurent property 
for a wide variety of recurrences, mostly (but not all) of the form (|2.8p . In particular, Fomin 
and Zelevinsky there gave the first proof of the Laurent property for the octahedron (discrete 
Hirota) recurrence (|2.7p . Subsequently, Speyer has developed a combinatorial model to prove 
more detailed properties of the Laurent polynomials generated by this recurrence - in particular, 
that all the coefficients are 1 [51] . 
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So far we have discussed the integrality of the sequence (j2.5p . but not the integrability of 
the Somos-4 recurrence. Taking {xQ,xi,X2,x-i) as coordinates, the map — > corresponding 
to ()2.4p preserves the degenerate Poisson bracket defined by 



{x„i,Xn}2 = {n-m)xmXn, (2.9) 
which has Casimirs 

un = ^7^^"+^ , n = l,2. (2.10) 

This bracket is of the 'log-canonical' type that has previously been found in the context of cluster 
algebras [18]; it is natural to consider it as a Poisson bracket on the field of rational functions 
C(xo, xi, X2, xs). (The reason for the subscript 2 on the bracket will become apparent in the 
next section.) The set of solutions of (12. 4p is invariant under the two-parameter Abelian group 
of gauge transformations generated by 

Xn^AXn, Xn^B'^Xn, A, B e C* . (2.11) 

The Hamiltonian vector fields corresponding to these transformations are respectively generated 
by the rational monomials 

K, = ^, K2 = (2.12) 

X2 X2 

which satisfy 

{Ki,Xn}2 = KiX„, {K2,Xn}2 = nK2Xn, { iiTi , K2 } 2 = A'l -f^2 . 



In fact, the most interesting part of the dynamics generated by (j2.4p takes place in the plane 
spanned by the Casimirs ui, U2 for the bracket {•, •}2- If we take the definition (|2.10p to hold 
for all n, then the quantities n„ are clearly invariant under the gauge transformations (|2.1ip . 
and satisfy the second-order recurrence 

aun+i + P , . 

Un+2 = W- (2-13) 

Un{Un+iy 

(So the fourth-order equation (j2.4p is the Hirota bilinearization of (j2.13p . which is a second-order 
ordinary difference equation.) By taking {ui,U2) as coordinates in C^, this corresponds to the 
rational map of the plane given by 

U2 

U2 ) V {au2 + (i) / {ui{u2f) 
which preserves the Poisson bracket 

{Un,Un+l} = UnUn+l, (2.14) 

or equivalently the symplectic form 

LOn = {UnUn+l)~^ du„ A dUn+1 (2.15) 

such that ujn+i = "^n- Furthermore, this has the conserved quantity 

/ 1 1 \ (3 , , 

J = Un Un+i + a \ \ H , (2-16) 

\Un Un+lJ UnUn+l 
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which defines a quartic curve 

CV-^Cr/ + a(C + ??) + /3 = 



(2.17) 



of genus one. Hence we see that ()2.13p produces a Liouville integrable system with one degree 
of freedom, and the curve (|2.17|) itself defines the two- valued correspondence u„ i— > Un±i, which 
is a particular case of the Euler-Chasles correspondence (see [57 1 l58 l [59]). 

Upon uniformizing the elliptic quartic we find that the explicit solution to ()2.13p is given by 



Un = p{z) - p{zo + nz), 
in terms of the Weierstrass p function for the elliptic curve 



E 



y2 =4X3-52X-<73, 



(2.18) 



(2.19) 



with 92 = 12A2 - 2 J, gs = AX^ - c/sA - a, A = {J'^/A - /?)/(3a), and zq, z G C/A = Jac(^) are 
given by elliptic integrals obtained from inversion of the relations p{z) = A, p(^o) = X — uq. The 
coefficients a, (3 and also J are given as elliptic functions of z by a = p'{z)'^, [5 = p'{zY {p{'^z) — 
p{z)),J = p"{z). 

From this it follows [27] that the solution to the initial value problem for the Somos-4 recur- 
rence ()2.4p can be written in terms of the Weierstrass sigma function as 



AB 



„ a {zo + nz) 
ai^zy 



(2.20) 



for suitable A, B. There is an analogous formula for the general solution of the Somos-5 recur- 
rence [28], which has an additional dependence on the parity of n. 
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Figure 1. A family of orbits for the nonlinear recurrence (|2.13p associated with Somos-4. 



The map defined by (j2.13p is a very simple example of the QRT family [l5]. It has a 2 x 2 
discrete Lax pair given by 



(2.21) 



where 



Ln('y) 



V - Un - Un+l 



-a 



a{l/Un + l/Un+l) + (3/ {UnU. 



ri+l) 
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M„ = , 

The equation (|2.13|) arises as the compatibility condition L„+iM„ = M„L„ for the system (|2.2ip . 
and the associated spectral curve is given by 

det (L„(i;) - vl) = v"^ - Ju + av + (3 = (2.22) 

where J is the conserved quantity given by (j2.16p . From the formulae v = ^r], v = Q + r] the 
elliptic quartic curve (|2.17|) is seen to be a ramified double cover of this rational (genus zero) 
spectral curve. 

It is clear from the above considerations that the dynamics of ()2.13p corresponds to a sequence 
of points P^ + nP on the elliptic curve E given by (j2.19p . or to the equivalent discrete linear flow 
zq + nz on its Jacobian, and in that sense (as was noted in [27]) it is the same as the underlying 
dynamics of the BT for the one-particle Garnier system constructed in [23] , or that of the BT 
for the g = 1 odd Mumford system as in [34]. However, while one can make changes of variables 
between (j2.13p and each of the latter two BTs, they are not canonical transformations, because 
the Poisson bracket (I2.14|) is incompatible with the Poisson structures of either of these BTs. 
Nevertheless, just as for the BTs, the recurrence ()2.13p is a discretization of a continuous time 
integrable system with the same Poisson structure and conserved quantities (in this case, only 
one of them), namely the flow in the plane with Hamiltonian J defined by p.l6p with n = 1, i.e. 

^ = {J,ui} = {aui + P)/u2 - {uifu2, 
at 

^ = {J,U2} = -{au2 + P)/ui + {U2fui. (2.23) 

From the same uniformization of the quartic (j2.17p as before, the solution of the system (I2.23P 
can be written down as 

ui{t) = p{z) - p{zo + p'{z)t), U2{t) = p{z) - p(zo + z + p'{z)t), 

and upon comparison with (j2.18p it can be seen directly how the discrete flow interpolates the 
continuous one (cf. Fig. [T]) . 

The construction of a sequence of points Pq + nP on elliptic curve E from a Somos-4 or 
Somos-5 sequence was previously understood in unpublished work of several number theorist^ - 
see the discussion of Zagier [67], and the results of Elkies quoted in [3]. The algebraic part of 
the construction is described in the thesis of Swart [53] (who also mentions unpublished results 
of Nelson Stephens), and van der Poorten has recently presented another construction based on 
the continued fraction expansion of the square root of a quartic [41 j . In fact, Somos-4 sequences 
have an ancestor from the 1940s, in Morgan Ward's work on elliptic divisibility sequences (EDS), 
which just correspond to multiples of a point nP € E [611 [62] i.e. this is the special case 
Pq = oo, so that zq = 0, with the further requirement that ^4 = i? = 1 in (I2.20|) . The iterates 
of an EDS, which are generated by (j2.4p with coefficients a = (^2)^, P = —X1X3 and integer 
initial data xi = 1, X2,a;3,X4 € Z with X2\x4, satisfy the divisibility property Xm\xn whenever 
m\n, and correspond to values of the division polynomials of the curve (for a description of 
these see Exercise 3.7 in [48]). In this sense, an EDS generalizes properties of certain linear 
recurrence sequences. For example, the Fibonacci numbers are generated by the recurrence 
Fn+i = Fn + Fn-i with initial values Fq = 1, Fi = 1, and form a divisibility sequence. Moreover, 
the even index terms Xn = i^2n form a divisibility sequence (so F2m\F2n whenever m\n) and also 
satisfy the Somos-4 recurrence 

F2n+A F2n-A = 9^2^+2 -F2n-2 — 8(^2^)^, 



For further references see http://www.math. wise. edu/'propp/somos. html 
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which corresponds to a degenerate case of the curve ()2.19p where the discriminant vanishes, so 
§2 — 27(7| = and the formula (|2.2U|) degenerates to an expression in terms of the hyperboUc 
sine. 

The arithmetical properties of EDS and Somos sequences - in particular the distribution 
of primes therein - are a subject of current interest [10\ \TT\ I49j . Some of these properties 
are discussed in the book [12] (see section 1.1.20, for instance), where it is suggested that 
such bilinear recurrences should be suitable generalizations of linear ones, with many analogous 
features. Based on the appearance of higher-order Somos recurrences in the work of Cantor 
on the analogues of division polynomials for hyperelliptic curves [6] (see also [ID] for analytic 
formulae), it was conjectured in [2^ that every Somos-A; sequence should correspond to a discrete 
linear flow on the Jacobian of such a curve (with an associated discrete integrable system), and 
the plausibility of this conjecture was justified by a nai've counting argument. However, on 
Propp's bilinear forum [H], Elkies had already given a more detailed argument to the contrary, 
based on a proposed theta function formula for the terms of such sequences, which indicated 
that while the general Somos-6 and Somos-7 sequence could be described by such a formula in 
genus two, the general Somos-A: for k > 8 could not. Thus in this setting the absence of the 
Laurent property appears to coincide with the absence of algebraic integrability. 

Nevertheless, in [1] it was shown that a particular family of solutions of Somos-8 recurrences 
can be described in terms of the Kleinian sigma function for a genus two curve (which is equiv- 
alent to an expression in theta functions), and these solutions are related to the BT for the 
Henon-Heiles system that was found in [231 [24] . The author has also found that the Somos-6 
and Somos-7 recurrences correspond to a rational map in with two independent conserved 
quantities, and there is a similar expression for the solutions in terms of genus two sigma func- 
tions. For instance, letting a denote the genus two Kleinian sigma function (see e.g. [5] for the 
definition), associated with a curve given by the affine equation = -|- c^x^ + • • • + cq with 
period lattice A, the expression 

x„ = AB"^<^51±5^. (2.24) 

where A,B G C*, vo,v £ mod A, satisfies a Somos-6 recurrence provided that v is con- 
strained according to 



1 1 1 

Pl2(v) Pl2(2v) Pl2(3v) 
P22(v) p22(2v) p22(3v) 



0, 



where Pjfc(v) = — 9^9^ log (t(v) are the corresponding Kleinian p functions. In the case of 
generic v, if this constraint does not hold, then Xn given by (j2.24p satisfies a Somos-8 recurrence 
instead. The full details will be presented elsewhere [3T] . 

Before moving onto other examples in the next section, we should mention one more feature 
of the Somos-4 recurrence, namely the fact that it generates solutions of a quartic Diophantine 
equation in four variables. If we rewrite the formula ()2.16p for the conserved quantity J in terms 
of the original variables Xn, we obtain the equation 

(a;„_i)^(3;„+2)^ + a (^Xn~-l{Xn+lf + {Xnfxn+2^ + /3 (XnfiXn+l)^ 

= J Xn-lXnXn+lXn+2- (2.25) 

If we have coefficients a,/3 G Z (or in Q), and if the Somos-4 recurrence (j2.4p with a set 
of integer initial data (xq, xi, X2, X3) generates a non-periodic sequence of iterates satisfying 
E Z for all n, then there are infinitely many quadruples of integers (x^-i, x„, x^+i, Xn+2) 
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that are solutions of the quartic Diophantine equation ()2.25p . (Note that in this CclSG, SiS long 
as all the integer initial data are non-zero, then the coefficient J which appears in (|2.25|) is 
uniquely determined, and J £ Q.) This can be seen as a particular instance of a general 
feature shared by all recurrences that both have the Laurent property and possess a rational 
invariant: generically, the orbit of suitable initial data will generate infinitely many solutions of 
an associated Diophantine equation. 

Diophantine Laurentness Lemma. Suppose that a kth-order rational recurrence of the 
form ()2.ip has coefficients in Q[c] (for some set of parameters c) and has the Laurent pro- 
perty, i.e. Xn G . . . ,x^^^,c] for all n. Suppose further that this recurrence also has 
a rational conserved quantity given by 

/l (^n ) ^n+l ) • • • ) 2^n+fe— 1 ) f') 

f2{Xm Xn+lj • • • ) — c) 

c] . // /2 / for some fixed integer values of c and initial 
data Xj = 1 or —1 for j = 0, . . . , k — 1, then the value of K G Q is fixed, and the recurrence 
generates infinitely many integer solutions of the Diophantine equation 

fl{Xn,Xn-\-l, . . . , c) = K /2(Xn, X^+l, • • • , c) 

as long as the corresponding orbit is not periodic. 

The integer sequence (j2.5p provides a concrete example of the above result: setting a = /3 = 1, 
the initial data 1, 1, 1, 1 yield the value J = 4 in (j2.25p . and for n > any four adjacent terms 
of this increasing sequence provide a distinct solution of the equation. In fact, in [S^ it is 
proved that the iterates of the Somos-4 recurrence satisfy the stronger property that Xn G 
Z[x^^ ,xi,X2,X3,a, (a^ + PJ)] for n > 0, which yields a broader set of sufficient criteria for 
integer sequences to be generated. In the next section we will see analogous results for some 
other recurrences. 



3 A fourth-order family 

A generalization of (j2.4p . that retains the Laurent property, is the family of fourth-order recur- 
rences 

Xn+2 Xn-2 = x'^+ix'^-l + x'^, (3-1) 

where the exponents a, b, c are positive integers. These generalized Somos-4 recurrences were 
first described in print by David Gale [T^, who noted that they all generate integer sequences 
from the initial values xq = xi = X2 = x^ = 1. Among various examples covered in [13], Fomin 
and Zelevinksy subsequently proved that each of these recurrences has the Laurent property. 
However, in contrast to the integrable structure of the original Somos-4 recurrence, most of these 
examples do not seem to correspond to completely integrable systems. 

Below we shall not present an analysis of the complete family (13. ip . but rather we focus on 
the special sub-family of recurrences defined by a = 6 = 1, with c G N. In this case, it will be 
convenient to introduce a parameter /3 as the coefficient of the second term on the right hand 
side; although this can always be removed by rescaling Xn, its inclusion preserves the Laurent 
property (while inserting another coefficient a in front of the Xn+sXn+i term does not, unless 
c = 2). These recurrences also satisfy the singularity confinement test that was proposed in [19] 
as an analogue of the Painleve test for discrete equations: if an apparent singularity is reached 
(in this case, corresponding to the situation that one of the iterates vanishes), then it is always 
possible to analytically continue through it. 
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Proposition 1. For each c G N the recurrence 

(3.2) 

which corresponds to the iteration of the rational map 



Xi X2 
I — > 

X2 Xs 

V X3 / \ {xiX3 + Px^)/xo ) 



(3.3) 



has the Laurent property in the sense that Xn S Z[x^^ , xf^ , x^^ , xf'^ , [3] for all n G Z, and also 
satisfies the singularity confinement test. Furthermore, (j3.3p is a Poisson map with respect to 
the log-canonical Poisson bracket {-j-jc defined by 

{xo,Xi}c = XqXi, {xo,X2}c = CXoX2, { Xq, X3} c = (c + 1)xoX3, 

{xi,X2}c = X1X2, {xi,X3}c = CX1X3, {x2,X3}c = X2X3, (3.4) 

which is nondegenerate for c ^ 2. 

Proof. The recurrence p.2p is of the cluster algebra type, so the Laurent property can be 
proved by the methods of [T^, where the details for the complete family (|3.ip are presented. 
However, here it is convenient to sketch a direct proof by induction, as this will have singularity 
confinement as an immediate corollary. The inductive hypothesis is that any four adjacent 
iterates Xk,Xk+i, Xk+2-, Xk+3 for < A; < n + 4 are coprime elements of the unique factorization 
domain TZ = Z[x^^, ^, x^^, x^^, /?], which has units ±Xq°, x^^ , Xg^Xg^ for £j G Z. Working 
mod x„+4, the congruences x„+5 = (3x'i^_^^x~l^, x„+6 = f3x''^_^^x~l2X~l^, x„+7 = (3"'^^ x"^'^ x'^^ 
all hold, so that 

Xn+&Xn+4: = ^n+3 ^xj^2^n+l^ {xn+3Xn+l + P ^n+2j = 0' 

since the bracketed expression in the middle is just x„+4X„ by (j3.2p . This proves the inductive 
step that Xn+8 £ and it is easy to see from (13. 2p that this element is coprime to Xn+5, 
Xn+6, Xn+7', the basc of the induction is trivial. This argument also demonstrates singularity 
confinement: if we have Xn+i = {xn+3Xn+i + /Jx^^j)/^™ = e ^ for some 72, so that x^^^g 
is potentially singular, then the preceding calculation shows that Xn+sXn+i = 0(e) and hence 
Xn+8 = 0{1) as e — > 0, so that the singularity is confined. The Poissonicity of the map ()3.3p 
is checked by a direct calculation, and in the coordinates y„ = logx„ the Poisson tensor for 
the bracket {^-jc is constant and has determinant (c — 2)^(c + 1)^. Thus for c G N it is 
nondegenerate unless c = 2, which gives the previously mentioned bracket (12. 9p preserved by 
the Somos-4 recurrence. ■ 



Remark 1. In the case of arbitrary parameters a = b and c, each of the recurrences (|3.2p admits 
a log-canonical Poisson bracket that generalizes {•, - jc; but there is no such bracket for a ^ b. 

Corollary 1. For each c the two- form 

' dxo A dxi dxo A dxg dx2 A dxg 
uj = I \ \ 

XqXi X0X3 X2X3 

dxoAdx2 dxiAdx3\ dxi A dx2 

H + (c + 1) , (3.5) 



X0X2 X1X3 / X1X2 



is preserved by the map (13. 3p . and this is symplectic for c 7^ 2. When c = 2 this two-form is 
degenerate, being the pullback of the two-form uJi in (j2.15p under the transformation 

(xo,Xi,X2,X3) 1-^ (mi,M2) = (xoX2/(xi)^, XiX3/(x2)^) . 
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We should now like to assert that the recurrences (j3.2p do not correspond to algebraically 
completely integrable systems when c > 3, based on the fact that in that case they have non- 
zero algebraic entropy. Recall that for a rational map, the algebraic entropy is defined as 
lim (log dn)/'n, where dn is the degree of the nth iterate of the map [21]. Usually finding this 

n— »oo 

limit requires extensive calculations of the corresponding sequence of rational functions of the 
initial data, or of the iterates of the projectivized form of the map. However, in this case we can 
exploit the fact that these recurrences generate Laurent polynomials, as well as the rescaling 
Xn P~^^^'^~'^^Xn (which for c 7^ 2 is equivalent to setting /3 = 1 in (j3.2p ). to argue that the 
degrees of the iterates as polynomials in the coefficient /? gives a suitable measure of the entropy. 

Proposition 2. For c G N, the nth iterate of (|3.2|) is a polynomial in jS of degree dn, as well 
as a Laurent polynomial in the initial data, where dn satisfies the recurrence 

dn+2 + dn-2 = max{ c d„ + 1, dn+i + dn-i } (3.6) 

for n > 2, with initial data do = di = d2 = d^ = 0. The algebraic entropy of the recurrence is 
zero for c = 0, 1, 2, while for c> 3 it is given by 

1 f c + \J(? - 4\ 

lim (log(i„)/n = - log 7, • (3-7) 

n-^oo Z \ Z / 

Remark 2. The full analysis of the 'tropical' (or piecewise- linear) recurrence (j3.6p is somewhat 
involved, and is omitted here, but we can mention that the determination of the value (j3.7p for 
the algebraic entropy follows from the fact that when c > 3 the degrees just satisfy the linear 
recurrence 



dn+2 + dn-2 = C(i„ + 1 

when n > 6, and hence they grow exponentially with n. For c = 0, 1 the growth of dn is linear 
in n, while for c = 2 it is quadratic in n (corresponding to the quadratic growth of logarithmic 
heights on elliptic curves [H]). Very similar analysis shows that for c > 3 the recurrences (13. ip 
fail Halburd's Diophantine integrability criterion [20], which requires that the logarithmic heights 
of all rational-valued iterates should grow no faster than polynomial in n. For instance, with 
initial data = x\ = X2 = xj, = \ each recurrence generates polynomials in and upon 

evaluating these at generic values of /3 G Q it can be demonstrated that the logarithmic heights 
of these numbers grow like the degrees dn- 

Having isolated the cases c = 0, 1, 2, we shall describe their integrable structure (in descending 
order). The case c = 2 is the original Somos-4 recurrence (|2.4p that was treated in the previous 
section, so we proceed with c = 1. 



Theorem 1. The map (]3.3p for c = 1 is superintegrable, in the sense that it has three indepen- 
dent conserved quantities Jk, k = 1,2,3, which satisfy 

{Ji,J2}i = = {Ji,J3}i, (3.8) 

where 

Ji = C0C1C2 — Cq — Cf — C2 + 2, J2 = Cq + Ci + C2-, = C0C1C2. (3.9) 

with 

XqX^ + x\ + X2 ^ XqxI + xfxs + XqxI + f3xiX2 

l-O — ! — , 

X1X2 XQX2X'i 
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^ X0X1X3 

The iterates of the corresponding recurrence 

Xn+2Xn-2 = Xn+lXn-l+ (iXn, (3-11) 

also satisfy the ninth-order linear recurrence 

X„+9 - (Ji + l){Xn+Q - Xn+s) - 3;„ = 0, (3.12) 

anc? the solution of the initial value problem for (13. lip has the form 

X3n+j = Aj TniJi/2) + Uni Ji/2) + /3Cj/(Ji - 2), j = 0, 1, 2, (3.13) 

where Tn and Un are the Chebyshev polynomials of the first and second kind respectively, and 
the coefficients Aj, Bj are given by 



_2x^_2pqiJ^-l^ 2x,+3 + /3C, . . 

Aj-2x, Ji(Ji-2) ' J, ^'^■'''> 



for j = 0,1,2. 

Proof. The proof of the above result is only sketched here, as further details will be presented 
elsewhere [31]. The main observation is that the recurrence (13. lip is linearizable, in the sense 
that the iterates satisfy the higher-order linear recurrence (j3.12p for a suitable Ji. (In the case 
P = 1, the nonlinear recurrence was originally considered by Dana Scott [17], who found that 
an integer sequence was generated from the initial data xq = xi = X2 = X3 = 1; in that case the 
linear recurrence (13.12P is satisfied with = 9.) In general one can take (|3.12p as the definition 
of iTi, and use (j3.1ip to back-substitute and rewrite it in terms of four adjacent iterates 

{Xn,Xn+l,Xn+2,Xn+3) = {p,q,r,s), 

as 

{p"^ + s'^)qr + (3{p + s){q'^ + r'^ +ps) + P'^qr 

Ji = , (3.15) 

pqrs 

which is found to be invariant with n, and defines a quartic threefold in C^. As a consequence, 
the linear recurrence (j3.12p holds for all n, and further implies the inhomogeneous linear equation 

Xn+6 - Jl Xn+3 + Xn + PCn = 0, (3.16) 

where the quantity C„ varies with n mod 3. Writing everything in terms of coordinates 

{xo,Xi,X2,X3) e & 

for the map (|3.3p with c = 1, this gives three independent quantities Cj given by (|3.1Up such 
that 

{Ji,Cj}i = 0, j = 0,1,2. 

These Cj are not preserved by the map, but symmetric functions of them are, which produces 
the formulae p.9p for three independent conserved quantities. The solution of the initial value 
problem can be conveniently expressed in the form (j3.13p . upon noting that (by separating 
out n mod 3) the homogeneous form of (j3.16p is equivalent to the second-order linear difference 
equation satisfied by the Chebyshev polynomials r„(^7i/2) = cos{n6), Un{J'i/2) = sm{n6)/ sinO 
with Jl = 2 cos 6*. ■ 
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Remark 3. The situation whereby an integrable system has more independent conserved quan- 
tities than the number of degrees of freedom is known as non-commutative integrability (in the 
sense of Nekhoroshev [37j), because not all these quantities can be in involution with one an- 
other. In this example, J\ Poisson commutes with both J2 and ^3, but {^j2iJz\\ 7^ 0. The 
terminology 'superintegrable' is applied in the even more special situation that the number of 
independent integrals is one less than the dimension of the phase space |64j . as is the case here. 



Upon applying the Diophantine Laurentness Lemma to the case of initial data (1, 1, 1, 1), and 
choosing integer /3 (with /? 7^ to avoid the degenerate case of a fixed point) we get infinitely 
many solutions of certain Diophantine equations corresponding to the conserved quantities. 

Corollary 2. With Ji = Ji{p,q,r,s) as in (j3T5]) . 

pr^{q + r + s) + sq^{p + q + r) + ps{pr + ps + qs) + Pqr{q + r) 

J2 — , 

pqrs 

{ps + + r'^){pr'^ + ps"^ + q^s + j3qr){pr'^ + q^s + j3qr) 



J3 



there are infinitely many integer solutions {p, q, r, s) of the double pencil of Diophantine equations 
given by 

Ai Ji + X2J2 + A3 J3 = Ai(/?2 + 6/? + 2) + A2(2/3 + 9) + SXaiP + 3)^, (3.17) 
for all f3£Z\ {0}, and any (Ai : A2 : A3) G 

As was remarked after Proposition El the Laurent polynomials generated by the c = case 
of (j3.2p show linear degree growth, so it might be anticipated that this case should also be 
linearizable. This indeed turns out to be so: the main results are very similar to the case c = 1, 
and are stated below without proof. 

Theorem 2. The map (j3.3p for c = is superintegrable, in the sense that it has three indepen- 
dent conserved quantities Jj, j = 1,2,3, which satisfy 

{JuJ2}o = = {Ji,J3}o, (3.18) 
and also { J'l, Qj }o = for j = 0,1,2, where 

Ji = Q0Q1Q2 -Qo-Qi- Q2, J2 = QoQi + Q1Q2 + Q2Q0 - 3, 

J3 = QoQiQ2, (3.19) 

with 

^ Xq + Xo ^ Xl + Xs ^ Xl Xo 

Qo = — Qi = - -, Q2 = — + — + ^^. 3.20 

Xl X2 Xq X3 X0X3 

The iterates of the corresponding recurrence 

Xn+2Xn-2 = Xn+lXn-l+ P, (3.21) 

satisfy the sixth- order linear recurrence 

Xn+6 - Jl Xn+3 + Xn = (3.22) 

and the solution of the initial value problem for (|3.21|) can be written explicitly in terms of 
Chebyshev polynomials of the first and second kind (Tn and Un respectively), as 

x^n+j = Aj Tn{Ji/2) + Un{Ji/2), j = 0, 1, 2, (3.23) 
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where the coefficients Aj, Bj are given by 

A- = 2x, - = -X, + ^ (3.24) 

for j = 0,1,2. 

We can apply the Diophantine Laurentness Lemma once more, taking initial data (1, 1, 1, 1) 
and /3 E Z, (with /3 / 0, — 1 to avoid periodic orbits), to get infinitely many solutions of quartic 
Diophantine equations corresponding to the conserved quantities for (I3.2ip . 

Corollary 3. With the identification {w,x,y,z) = (x„, x^+i, Xn+2, a^n+s), when n = the 
relations (j3.19p define Jk = Jk{w, x, y, z) for k = 1,2, 3, as 

~ {w"^ + z^)xy + [5{wx + wz + yz) 
wxyz 

+wy{x'^ + y'^ + z^ + wy) + I3{x'^ + y'^ + xz + wy) ) '^'^^y^>' 

~ (u; + y)(x + 2;)(x2; + yif; + /3) 

Js — 1 

wxyz 

and there are infinitely many integer solutions {w, x, y, z) of the double pencil of Diophantine 
equations given by 

Ai Ji + \2J2 + As = Ai(3/? + 2) + A2(4/3 + 9) + 4A3(/3 + 2), (3.25) 
for all p eZ\ {0, -1}, and any (Ai : A2 : A3) G P^. 

Remark 4. The initial data (1,1,1,1), together with the restrictions on /3, are sufficient to 
ensure that the each of recurrences (j3.1ip and (j3.2ip generate non-periodic integer sequences, and 
hence infinitely many solutions of the corresponding Diophantine equations, given in Corollary [2] 
and Corollary [3] respectively. However, due to the fact that the recurrences are integrable (and 
even linearizable) in both cases, it is possible to choose much more general initial data and still 
generate integer sequences, which produce solutions of the same Diophantine equations but with 
different values on the right hand sides of (|3.17p and (|3.25p respectively. 



4 Outlook 

The Laurent property appears to be an extremely elegant, but somewhat special, feature of 
certain rational maps. In particular, it seems to hold for integrable bilinear or discrete Hirota 
type equations, such as (j2.4p and (j2.7p . but also for the whole family (j3.2p . whose members have 
non-zero algebraic entropy for c > 3. For the latter family, we have noted the close connection 
between the Laurent property and the notion of singularity confinement as introduced in |19j . 
(For other examples of confined maps with the Laurent property see [291 130j.) This connection 
seems to persist for rational maps that do not themselves have the Laurent property. 
For example, consider the second-order equatiord 

K)' + 1 

Un+l = , (4.1) 

which is superficially very similar to (j2.13p . and preserves the same symplectic form (j2.15p . 
The real phase portrait in also looks qualitatively similar: Fig. [2] seems to display the same 



am grateful to Vasilis Papageorgiou for showing me this example, which I believe is due to Claude Viallet. 
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structure of invariant curves as Fig. [TJ Furthermore, the equation (|4.ip satisfies singularity 
confinement, with the singularity pattern being e, e~^, e~^, e~^, e (for e — > 0), which suggests 
making the substitution 



Tn+2Tn-2 



T„+i(T„)2r„_i ■ 

Thus Un given as above satisfies ()4.ip whenever r„ satisfies 

Tn+3 7"n-3 = ('^«+2Tn-2)^ + (^n+l )^ (t„)'^(t„_i )^ , 



(4.2) 



(4.3) 



and this sixth-order recurrence has the Laurent property, as well as satisfying the singularity 
confinement test. The singularity pattern for (|4.ip . which includes poles, "unfolds" to yield 
isolated zeros, i.e. r„ = e for some n with adjacent iterates being 0(1) as e — > 0. 

However, the logarithmic heights h{Tn) of rational iterates grow exponentially with n. To 
see this, it is instructive to take all six initial values for (j4.3p equal to 1, yielding the integer 
sequence 

1, 1, 1, 1, 1, 1, 2, 5, 29, 1241, 3642581, 80305336110269, . . . , 
which grows like log /i(t„) = loglog|T„| ~ nlog7 with log 7 ~ 0.733, where 7 ~ 2.081 is the 



largest modulus root of the polynomial 7 



7 



27 —7 + 1. Hence the logarithmic heights 



of the rational numbers ti„ that lie in the orbit of {uq,ui) = (1,1) generated by (14. 1|) grow 
exponentially, and Halburd's Diophantine integrability criterion is failed. Similar arguments 
hold for generic orbits, and it follows that the curves appearing in Fig. [2] are not algebraic. To 
see this, recall that by the Hurwitz theorem a curve with an infinite order automorphism group 
has genus zero or one, and under iteration of such automorphisms the logarithmic heights of 
rational points grow linearly on a rational curve and quadratically on an elliptic one 




Figure 2. A family of orbits for the nonlinear recurrence (|4.1 



We have concentrated on recurrences of the particular form (j2.3p . but this is not necessary 
for the Laurent property. Another interesting (and algebraically non-integrable) example is the 
second-order equation 



Un+l + Un-l = Un + 



(n„,)2' 



a/0. 



(4.4) 
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which in [21] was found by Hietarinta and Viahet to display singularity confinement with the 
pattern e, e~^, e~^, e, yet it has positive algebraic entropy and its real orbits display the char- 
acteristics of chaos. By way of the substitution 

_ Tn+2Tn-l 
Un — / NO 
{Tn+lTnY 

we arrive at the fifth-order recurrence 

Tn+z{Tnf{Tn-lf = {T^+^f {Tn-lf - {Tn+2? {Xn+lf Tn-2 + a(r„+l-^„)^ (4.6) 

which itself satisfies singularity confinement and has the Laurent property, i.e. for all n the 
iterates satisfy Tn G I^It^^ , . . . ,r^^,a]. In this case, the logarithmic heights of rational iterates 
r„ G Q (with a G Q) generically satisfy /i(r„) ~ CC" for some C > and C = (3 + a/5)/2 is 
the square of the golden mean, while logC turns out to be the value of the algebraic entropy 
for (j4.4p . Note that while the calculation of the algebraic entropy is quite involved |2H [55] , it 
is quite straightforward to calculate the growth of heights from ()4.6p . Similarly to the previous 
example, the only confined singularities that appear in (j4.6p are isolated zeros. 

These examples illustrate the following general phenomenon: whenever we have a rational 
map with confined singularities, including poles, it should always be possible to "unfold" these 
into confined zeros, by embedding the map in higher dimensions via a change of variables, 
and the new map thus obtained should have the Laurent property. This is analogous to the 
way in which continuous integrable systems with the Painleve property, that have meromorphic 
solutions, admit a Hirota bilinear form (or multilinear form) in terms of tau- functions that are 
holomorphic. Although this phenomenon (the existence of a tau-function) is very well known 
for discrete integrable systems [46], so far it does not seem to have been to have been exploited 
in the case of non-integrable maps. 



(4.5) 
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